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ABSTRACT 

We apply the disc instability model for soft X-ray transients to identify system pa- 
rameters along evolutionary sequences of black hole X-ray binaries (BHXBs) that are 
consistent with transient behaviour. In particular, we focus on the hitherto neglected 
group of BHXBs with intermediate-mass giant donor stars. These spend a significant 
fraction of their X-ray active phase crossing the Hertzsprung gap. 

Three Case B binary sequences with a black hole accretor and 2.5M Q initial donor 
mass are presented in detail. We formulate rules which summarize the behaviour of 
these sequences and provide an approximate description for Case B mass transfer in 
intermediate-mass BHXBs. Chiefly, the timescale of the overall radius expansion is 
given by the initial donor mass, while the surface appearance is determined by the 
current donor mass. 

With these rules we obtain a general overview of transient and persistent behaviour 
of all intermediate-mass BHXBs by just considering single star sequences of different 
mass. We find that although systems in the process of crossing the Hertzsprung gap 
are in general persistently bright, with Eddington- or super-Eddington transfer rates, 
there is a narrow instability strip where transient behaviour is possible. This strip 
extends over a secondary mass range 2.0 Ss M% Ss 3.5M©. GRO J1655-40 might 
be such a system. We predict that there are no BHXB transients with donors more 
massive than 3.5Mq, and no neutron star transients in the Hertzsprung gap. 

Key words: accretion, accretion discs — binaries: close — black hole physics - 
stars: evolution — stars: individual (GRO J1655-40) 



1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) are semi-detached com- 
pact binaries with a black hole (BH) or neutron star (NS) 
primary that accretes mass from an accretion disc. A Roche 
lobe filling low-mass main-sequence or (sub)giant star feeds 
the disc. LMXBs appear in two main varieties: They are ei- 
ther persistently bright, or have outbursts and long phases 
of quiescence (e.g. White & van Paradijs 1996). The lat- 
ter, known as soft X-ray transients (SXTs; cf. Tanaka & 
Shibazaki 1996 and Chen et al. 1997 for reviews), are be- 
lieved to have an unstable accretion disc which alternates 
between a hot and cool state (Cannizzo et al. 1982, Taam 
& Lin 1984, Cannizzo, Chen & Livio 1995, King & Ritter 
1998), similar to dwarf novae among cataclysmic variables. 

A simple criterion for stable disc accretion is that the 
temperature Td at the outer disc edge is larger than the 
hydrogen ionization temperature Th — 6500 K (e.g. King, 
Kolb & Burderi 1996; hereafter KKB). As Td increases with 
the mass transfer rate M, a system is persistently bright if 



M is larger than the critical rate M cr it where Td = Th- The 
system is transient if M < M cr i t . 

An important difference between dwarf novae and SXTs 
is the dominant role irradiation of the accretion disc by 
the central accreting source plays in LMXBs. Van Paradijs 
(1996) showed that the observed critical transfer rate sep- 
arating transient from persistent neutron star LMXBs is 
much lower than the critical rate separating dwarf novae 
from persistently bright (novalike) cataclysmic variables. 
This can be understood quantitatively by assuming that ir- 
radiation dominates over viscous heating at the outer disc 
rim, hence stabilizing the disc even at a smaller mass transfer 
rate. Similarly, irradiation naturally explains the observed 
large ratio of optical to X-ray flux from the disc and the 
observed long, quasi-exponential decline after an SXT out- 
burst (King & Ritter 1998). 

Models for the evolution of LMXBs predict the mass 
transfer rate and thus by comparison with M cr it the ap- 
pearance of the system as a transient or persistently bright 
X-ray source (assuming that the instantaneous transfer rate 
is close to the evolutionary mean transfer rate). This is a 
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powerful diagnostic tool for testing evolutionary models, or, 
in turn, the disc instability model for SXTs. King, Kolb & 
Szuszkiewicz (1997; hereafter KKS), King et al. (1997) and 
KKB investigated short-period systems with essentially un- 
evolved donors, and long-period systems having low-mass 
giant donors with degenerate helium cores. In the following 
we extend this work, and present evolutionary models for 
black hole X-ray binaries (BHXBs) with intermediate-mass 
donors in the process of crossing the Hertzsprung gap, a 
hitherto largely neglected parameter range. We explore to 
what extent such systems appear as SXTs. A major moti- 
vation for this study is provided by recent optical observa- 
tions of the transient X-ray source GRO J1655-40 (Orosz & 
Bailyn 1997), well-known from the apparent superluminal 
motion in its radio jets (Hjellming & Rupen 1995). System 
parameters derived for GRO J1655-40 place this BH binary 
just in the above parameter range dicussed here. Some of 
the results presented below have already been applied to a 
discussion of the evolutionary state of GRO J1655-40 (Kolb 
et al. 1997). 

We put intermediate-mass BHXBs into context by re- 
viewing main concepts of LMXB evolution briefly in Sect. 2. 
In Sect 3 we present in detail three evolutionary sequences, 
with parameters close to the ones observed in GRO J 1655- 
40. We generalize the results of these calculations in Sect. 4 
in the form of three rules which approximately describe the 
general evolution of intermediate-mass BHBXs. These rules 
allow one to obtain a general overview of transient and per- 
sistent behaviour of all intermediate-mass BHXBs by just 
considering single star sequences of different mass. Sect. 5 
concludes with a critical discussion. 



2 EVOLUTION OF LMXBS 

The mechanism driving mass transfer naturally divides 
LMXBs into two distinct classes. In the first group, nuclear 
expansion of the secondary maintains the semi-detached 
state, while in the second group, orbital angular momentum 
losses such as gravitational radiation and magnetic braking 
drive mass transfer. For convenience we denote neutron star 
or black hole LMXBs in the two groups with n-driven or 
j-driven systems. 

Generally, n-driven LMXBs are long-period systems 
with (sub)giant donor stars on the first giant branch. Hy- 
drogen burns in a shell source above a nuclearly inactive 
helium core. The orbital period increases until mass trans- 
fer terminates when core helium burning ignites, or when 
the hydrogen-rich envelope is fully transferred to the pri- 
mary, whichever happens earlier. The final state is a wide, 
detached binary with a white dwarf secondary and a BH 
or NS primary (possibly a millisecond pulsar). The LMXB 
lifetime depends on the initial secondary mass and orbital 
distance and can be as long as 10 9 yr if both are small, but 
is typically 10 7 - 10 s yr. 

Conversely, j-driven systems are short-period binaries 
with main-sequence donors and evolve to shorter orbital pe- 
riod. If the j-driven evolution is the same as for cataclysmic 
variables (see e.g. Kolb 1996, King 1988 for reviews) the 
bright phase where magnetic braking operates (P i; 3 h) 
typically lasts 10 s yr. For P <, 3 h mass transfer continues 
indefinitely at a lower rate, but the period increases again 



when the secondary becomes degenerate, i.e. when its mass 
is <, 0.06M Q . 

A critical bifurcation period Pt, ~ 1 — 2 days (Pyl- 
yser & Savonije 1988, 1989) separates j-driven and n-driven 
LMXBs. Systems born with P < Pb evolve towards shorter 
periods, and those born with P > P b towards longer P. Sim- 
ilarly, the initial donor mass Mi\ plays a role in determining 
the group membership: If Mi\ Ss 0.8Mq the system is j- 
driven, as the donor's main-sequence lifetime is longer than 
a Hubble time. If M-z\ ii 1.5Mq the system cannot evolve to 
shorter orbital periods as magnetic braking does not oper- 
ate in such massive stars, and gravitational radiation alone 
is too weak to dominate the evolution. In the intermediate 
range 0.8 & 1.5 both groups exist. Mass transfer sta- 

bility defines an upper limit for the mass ratio q = M2/M1 
(M2 is the donor mass, Mi the primary mass), hence for 
Mi and M^. This limit depends on the donor's response to 
mass loss (Hjellming 1989), i.e. is a non-trivial function of its 
evolutionary state. As a rule, the limit is roughly M2 SS M\ 
in the case of conservative mass transfer, but can be some- 
what larger with mass loss from the system (cf. Kalogera 
& Webbink 1996, where the NS case is discussed in detail). 
If the donor is more massive the system suffers a short, vi- 
olent phase of mass transfer/loss and is unlikely to appear 
as an X-ray source. Obviously, neutron star X-ray binaries 
with Roche-lobe filling donors more massive than 2 — 3M Q 
do not exist, while there is no such upper limit for BHXBs, 
hence no clear separation between low-mass and high-mass 
BHXBs. The importance of BHXBs with intermediate-mass 
donors was first pointed out by Romani (1994). 

Recently, two limiting cases of LMXB evolution have 
been considered to determine the incidence of transient be- 
haviour among LMXBs: KKB and KKS described the j— 
driven evolution of completely unevolved (ZAMS) donors. 
They found that BH systems are always transient, consis- 
tent with observations, and that donor stars in NS systems 
need to be very close to the end of core H burning in order to 
be transient, consistent with the observed small fraction of 
transient NS LMXBs. King et al. (1997) and KKS studied n- 
driven LMXBs with low-mass donor stars well-established 
on the first giant branch, and found that essentially all of 
them are transient. They considered giant donors with a thin 
H burning shell source above a degenerate He core where 
the stellar radius R and luminosity L are unique functions, 
R oc M c 6,1 , L oc M 8 ' 1 , of the slowly growing core mass M c , 
essentially independent of the total stellar mass (e.g. Web- 
bink, Rappaport & Savonije 1983). Stars with mass S5 2Mq 
indeed establish such a structure soon after the end of cen- 
tral H burning and spend a considerable time ascending the 
first giant branch, before the ignition of helium burning in 
the centre terminates their radius expansion. Most n-driven 
NS LMXBs are well described in this way, but not BHXBs 
with more massive donor stars (M <; 2.0 — 2.5M Q ). These 
ignite core helium burning before the core becomes highly 
degenerate, i.e. before the first giant branch evolution gov- 
erned by the above core mass relations takes hold. Simplified 
descriptions for the evolution of these stars fail, and full stel- 
lar models are needed. 

We consider these n-driven intermediate-mass BHXBs 
in the next sections. 
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3 ILLUSTRATIVE CASE B MASS TRANSFER 
SEQUENCES WITH BLACK HOLE 
ACCRETORS 

The binary evolutionary phase we consider here is Case B 
mass transfer: the donor star fills its Roche lobe after termi- 
nation of core H burning but before the ignition of core He 
burning. Case B mass transfer is well-studied (e.g. Kippen- 
hahn et al. 1967, Pylyser & Savonije 1988, Kolb & Ritter 
1990, De Greve 1993) in the context of Algol evolution. All 
these studies assume that the donor is more massive than 
the accretor, simply because the faster evolving more mas- 
sive component always fills its Roche lobe first. In BHXBs, 
with a more complicated evolutionary history, Case B mass 
transfer begins with an already inverted mass ratio, i.e. the 
donor is less massive than the BH. Hence there is no initial 
rapid phase with extremely high transfer rate. 

In the following we present calculations of three Case 
B mass transfer sequences on to a BH primary from a sec- 
ondary with initial mass close to the transition region be- 
tween degenerate and non-degenerate helium ignition. We 
report briefly on the numerical technique and associated 
problems (Sect. 3.1), describe in detail the evolution of the 
secondary star (Sect. 3.2), and then focus on the resulting 
appearance in X-rays along the sequences in Sect. 3.3. 



3.1 Input Parameters and Computational 
Technique 

We used Mazzitelli's stellar evolution code in a version as de- 
scribed by Mazzitelli (1989; see also references therein) with 
pre-OPAL opacities. The computations started from chemi- 
cally homogenous ZAMS models with Population I mixture 
(X — 0.70, Y — 0.28). Convection is treated by the standard 
mixing length theory; a calibration to a solar model deter- 
mines the mixing length parameter to 1.4 (no overshooting 
was allowed). 

To allow the application to binary evolution the mass 
transfer rate M was calculated for each time step according 
to 

^ K exp (m) (1) 

(cf. Kolb & Ritter 1990), where AR = R - R h is the dif- 
ference between the secondary's radius R and Roche lobe 
radius i?L, H the atmospheric pressure scale height, and k 
a numerical constant, usually taken as k = 1. Although (|l|) 
is strictly valid only for AR < and has to be replaced by 
an expression with an approximate power-law dependence 
on AR/H for AR > 0, we use (Jl]) for any AR. Furthermore, 
to avoid numerical instabilities at high transfer rates we set 
k = 10 for M ^ 1O" 7 M yr" 1 . This procedure is justified 
as long as AR/R < 1 (Kolb & Ritter 1990), i.e. as long as 
the stellar radius is close to the critical Roche radius. For 
the sequences presented here this is always the case. A high- 
resolution reference sequence with k = 1 was calculated to 
check the validity of sequences obtained with k = 10 ex- 
plicitly in the case of sequence S3 (see below). Attempts to 
retain the damping factor k = 10 also in later phases (where 
AR < 0) failed as H/R becomes non-negligible for extended 
giants. A damped evolution would proceed qualitatively dif- 
ferently from the true evolution with k = 1, leading to higher 
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Figure 1. Radius R (upper panel) and effective temperature T e fj 
(lower panel) of a 2.5Mq single star (dashed) and the secondary 
star in the binary sequences SI (dash— dotted) , S2 (thick solid 
line) and S3 (thin solid line), cf. Tab. 1, as a function of time 
elapsed since the (donor) star left the main sequence. 



Table 1. Binary sequence parameters: BH mass Mi, donor mass 
M2, orbital separation a; the last columns lists an internal se- 
quence number. 



sequence 


Mi/M Q 


M 2 /M Q 


a/R.0 


int 


SI (initial) 


6.8 


2.5 


13.2 


b860 


S2 (initial) 


6.8 


2.5 


15.2 


b849 


S3 (initial) 


6.8 


2.5 


18.5 


b830/32 


S3 (final) 


8.726 


0.574 


213.1 





mass transfer rates and a correpondingly earlier termination 
of mass transfer. 

3.2 Detailed Description of the Sequences 

All three binary sequences, S1-S3, begin mass transfer with 
a 2.5M0 donor star and a 6.8M0 black hole. The initial or- 
bital separation is a; = 13.2 (SI), 15.2 (S2) and 18.5Rq (S3), 
see the summary in Tab. 1. Mass transfer is conservative, 
i.e. the total binary mass and orbital angular momentum 
is constant. Sequence S3 follows the donor star until igni- 
tion of central helium burning. The subsequent evolution 
is detached, and we evolved the secondary further through 
central helium burning and the brief subsequent asymptotic 
giant branch phase with weak thermal pulses. In contrast, 
sequences SI and S2 have been terminated once the sec- 
ondary was established on the first giant branch. 

In Figs. 1 - 5 we show characteristic parameters along 
the sequences S1-S3, sometimes together with the evolu- 
tion of a 2.5M0 single star (for convenience referred to as 
sequence SO). 

Figure 1 (upper panel) reveals that the overall radius 
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Figure 2. Surface luminosity L (top), nuclear luminosity L nuc 
(middle panel) and absolute value of the gravo-thermal luminosity 
L g = L — L nuc (bottom) of a 2.5Mq single star (dashed) and the 
secondary star in binary sequences S1-S3, as a function of time 
elapsed since the (donor) star left the main sequence. The linestyle 
is as in Fig. 1. 



evolution with time is very similar in all sequences S0-S3, ex- 
cept for the "plateau phase" soon after mass transfer turn- 
on, which corresponds to a brief radius contraction phase at 
the red end of the Hertzsprung gap of SO. This phase occurs 
at a characteristic radius, the "plateau radius" R p , and will 
be examined in more detail below. The similarity shows that 
significant mass loss from the outer envelope has little effect 
on the overall expansion of the donor star, which is gov- 
erned by the nuclear evolution of the core region. The overall 
temporal evolution of the surface luminosity L is systemati- 
cally different for sequences with and without mass transfer 
(Fig. 2): in S1-S3 mass loss causes a steep build-up of gravo- 
thermal luminosity L g (L g < 0) and a slight reduction of 
the nuclear luminosity L nuc from the hydrogen shell source. 
This behaviour is superimposed on the luminosity evolu- 
tion pattern of SO in the Hertzsprung gap, where hydrogen 
burning changes from a thick to a thin shell. The net effect 
is that thin hydrogen shell burning is established earlier, at 
a slightly smaller core mass. Once thin H burning is fully 
established the luminosity grows with further growing core 
mass. The corresponding effective core mass-luminosity re- 
lation is much steeper than the standard relation L oc M,?' 1 
for fully degenerate cores, and meets the latter one only im- 
mediately before core helium ignition (Fig. 3). As a result, 
the secondary in S1-S3 begins its ascent along the first giant 
branch at a significantly smaller surface luminosity than the 
2.5Mq single star. The drop of L, a familiar feature of Case 
B mass transfer, goes along with a correspondingly fast de- 
crease of Teg to logr o fr/K ~ 3.70, i.e. the secondary in the 
mass loss sequence crosses the Hertzsprung gap even faster 
than the 2.5M Q single star (Fig. 1, lower panel). Roughly, 
at any time along S1-S3 in the Hertzsprung gap, L and T e g 



Figure 3. Surface luminosity L as a function of core mass M c 
(defined as the mass enclosed by the shell where the H burning 
energy generation rate is maximal). Full triangles: selected mod- 
els along sequence S3. The triangles with smallest and largest 
M c correspond to models at the luminosity minimum and close 
to central He ignition, respectively. The open triangle marks a 
2.5Mq single star at its luminosity minimum. The positions of 
models with larger M c along the single star sequence coincide 
with the positions of sequence S3. The solid curve is the core- 
mass luminosity relation for low-mass (< 1.5Mq) giant stars, cf. 
King ct al. 1997. 

are the same as for a single star with the same mass and 
radius as the secondary star at that time. 

With the ignition of helium burning the secondary's ra- 
dius decreases sharply and mass transfer terminates. At this 
point the binary parameters of S3 are Mi = 8.726M0, Mi = 
0.574M Q (with core mass O.32OM ) and P = 118.3 days, 
a = 213.1R.Q (see Tab. 1). For the next 270 Myr the sec- 
ondary burns helium in the centre and is well inside its 
Roche lobe (R ~ 8 - 10R© versus R L = 34.6R ). Subse- 
quently, He burning moves to a shell source and the radius 
increases again, at a core mass of ~ 0.420M Q . The ensu- 
ing double shell burning is unstable and leads to thermal 
pulses. We find that these pulses are rather weak, with de- 
creasing maximum radius. The radius maximum in subse- 
quent pulses is always smaller than the Roche radius (though 
not very much), and quickly decreasing. Eventually the sec- 
ondary contracts towards the white dwarf stage. If there are 
no significant wind losses during the thermal pulses the fi- 
nal state of S3 is a wide, detached binary with a ~ 9Mq 
black hole primary, a ~ 0.55Mq carbon-oxygen white dwarf 
companion and ~ 5.6 months orbital period. 

3.3 Transient Versus Persistent Accretion 

We distinguish transient from persistent BHXBs by compar- 
ing the mass transfer rate M with the critical rate M cr i t for 
which hydrogen just ionizes at the outer disc edge. A system 
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Figure 4. Mass transfer rate M as a function of time (elapsed 
since the star left the main sequence) for sequences SI (bottom 
panel), S2 (middle panel) and S3 (upper panel). The linestyle is 
broken when the mass transfer rate is smaller than the critical 
rate M cr i t for disc instability (taken from KKS). 



is transient if M < M cr n. In Fig. 4 we plot M as a func- 
tion of time along sequences S1-S3. The linestyle is broken 
if M < Merit, with Merit taken from KKS. Obviously, the 
accretion disc is stable for about ~ 4 — 5 Myr when the sec- 
ondary crosses the Hertzsprung gap, and the transfer rate is 
close to (or slightly above) the Eddington limit 

Mdd 2 x lO^Moyr" 1 (jj^J (2) 

(assuming an effective accretion efficiency rj = L acc /Mc 2 = 
0.1). The system becomes a transient source soon after it 
begins to ascend along the Hayashi line, and remains tran- 
sient until helium ignites in the centre ~ 12 Myr later. 
Significantly, there is also a short transient phase (lasting 
~ 1 Myr) right in the Hertzsprung gap, centered around 
log T e ff /K = 3.73. This phase is caused by the transition 
from thick to thin shell source burning, the same effect that 
leads to a brief radius contraction when the 2.5M single 
star crosses the Hertzsprung gap. The associated change of 
the internal structure of the star is highly non-linear and 
complex. Just as the ultimate cause for the expansion to 
the red giant state cannot be easily derived from a simple 
analysis of the stellar structure equations, there is no easy 
way to understand this temporary halt in the expansion ei- 
ther. 

The Hertzsprung gap transient phase in S1-S3 and the 
radius contraction phase of SO all occur at the same range 
of T e ff (Fig. 5), while the corresponding plateau radius R p 
is largest for SO, and progressively smaller the earlier mass 
transfer began in the Hertzsprung gap. The secondary mass 
in the plateau phase is 1.9, 2.05 and 2.25M© for sequences 
SI, S2 and S3, respectively. In fact, R P (M2) is practically 



Figure 5. Mass transfer rate M as a function of effective tem- 
perature T c ff for sequences SI (dash-dotted), S2 (thick solid line) 
and S3 (thin solid line). Also shown is the critical rate for disc 
instabilities (from KKS) for sequence SI. 

the same as the plateau radius of a single star with the 
secondary's mass in the plateau phase. 

Details of the exact location (e.g. the radius R p ), du- 
ration and depth of the radius contraction phase of a single 
star must necessarily depend on details of the stellar input 
physics, in particular on opacities and treatment of convec- 
tion. As an example, the 2.5M single star track obtained 
by Salaris et al. (1997) with OPAL opacities shows a much 
less pronounced radius contraction phase at a slightly higher 
effective temperature (logT e fr/K ~ 3.75) than our sequence 
SO. Hence we do not expect our models to reproduce the pre- 
cise location of the real transient phase in the Hertzsprung 
gap, and the actual mass transfer rate in this phase. But 
our sequences certainly show the differential change of that 
phase with varying initial orbital distance: the earlier mass 
transfer starts in the Hertzsprung gap, the less pronounced 
the decrease of M in the transient phase. If the secondary is 
already very close to the plateau phase when it fills its Roche 
lobe for the first time, the system might almost detach in the 
transient phase and not appear as an X-ray source at all; but 
if the system begins mass transfer early in the Hertsprung 
gap it certainly will. 



4 GENERALIZATION: 

INTERMEDIATE-MASS BHXB EVOLUTION 

We formulate three rules which summarize the behaviour of 
the above evolutionary sequences and provide an approxi- 
mate description for the Case B evolution of intermediate- 
mass BHXBs. 

(Rl) the overall time evolution R(t) of the donor star radius in 
a binary sequence with initial donor mass M2i is the same as the 
time evolution of the radius of a single star with mass Mo = M2i 
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while 

(R2) the plateau phase (the brief contraction phase or slow- 
down of the expansion phase) in the Hertzsprung gap occurs at a 
donor radius R p equal to the plateau radius of a single star with 
the donor's mass 



(R3) a donor star with mass M 2 and radius R has the same 
effective temperature T e g (and surface luminosity L) as a single 
star with mass Mo = M 2 and radius R 

In other words: the timescale of the overall radius ex- 
pansion is given by the initial donor mass, while the surface 
appearance and the plateau phase is given by the current 
donor mass. 

In the following we use these rules to investigate 
the incidence of transient and persistent behaviour in 
intermediate-mass BHXBs with arbitrary initial mass and 
initial separation, simply by considering the evolution of sin- 
gle stars in this mass range. 

In order to do so we note that in the case of stationary 
mass transfer (R/R = Pl/Pl) the transfer rate is given by 



M = Ah 



K 



Coq - CH 



(3) 



(cf. e.g. Ritter 1996). Here K = (dlnP/dt.) A /=const. is the 
radius expansion in the absence of mass loss (e.g. due to 
nuclear evolution), £ eq the thermal equilibrium mass-radius 
exponent, and £r the Roche lobe index (£r ~ 2M2/M1 — 5/3 
for conservative evolution). We set £ oq = and, using (Rl), 
approximate K by dlnP/di along the single star sequence. 
It is clear that these secondaries are not in thermal equi- 
librium (Fig. 2), and their radius expansion is not purely 
nuclear. However, the mass loss perturbs the overall radius 
evolution only slightly. 

Then we obtain the approximate binary evolution for 
given initial donor mass M 2 i by integrating M(t), i.e. K(t) 
of a single star with mass Mo — M<n, over time. To decide 
the stability of the accretion disc we need to know the mass 
transfer rate M as a function of orbital period P. The main 
characteristics of the real evolutionary tracks M(P) are al- 
ready evident from Fig. 6 where we plot the mass trans- 
fer rate Mi soon after turn-on of mass transfer, estimated 
from (^|), as a function of initial orbital period Pi, for sec- 
ondaries in the mass range of interest (M 2 i = 1.3 — 3.5M0). 
Roche geometry determines the initial period to log(Pi/d) = 
1.5 log(P/R ) - 0.5 log(M 2i /M ) - 0.433. Also shown is the 
critical rate Merit for disc instability, taken from KKS, and 
the Eddington accretion rate (Q) for a typical 8M0 BH ac- 
cretor. 

For M 2 i ^ 1.9M we find that Mi < Merit at any 
Pi, while for M 2 i <; 1.9M the transfer rate is smaller 
than M cr it either if the donor is in the plateau phase, 
P = P(P P (M2i),M 2 i) = P p (M 2i ), or if P is sufficiently 
long, i.e. longer than a critical period P cr it where Mi = 
min(M cr it, M c dd). The Eddington rate enters because for 
M cr it > Medd the mass accretion rate in any system with 
super-Eddington mass transfer rate is probably limited to a 
value near M e dd, hence the accretion disc is unstable even if 
M > Merit, (assuming that the character of disc irradiation 
does not change for super-Eddington mass transfer; note 
that because Merit is determined by irradiation, it is the 
central accretion rate which enters the stability criterion.) 



of 
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Figure 6. Transfer rate Mi versus orbital period Pi at turn- 
on of mass transfer, for three different donor masses (full line: 
3.5M ; dashed: 2.5M ; dotted: 1.3M ), estimated from Eq. (^) 
and assuming a 8M primary. The horizontal line indicates the 
corresponding Eddington transfer rate, see (J2J) - The thick solid 
line is the critical mass accretion rate M cr it separating transient 
(M < Mcr) from persistent (M > M cr ) systems, taken from KKS. 



We show Pcrit as a function of M 2 i in a plot of orbital 
period versus donor mass (Fig. 7; a similar figure was al- 
ready published in Kolb et al. 1997). This critical period 
separates BHXBs with stable accretion disc (unhatched re- 
gion) from systems with unstable discs (hatched region) at 
turn-on of mass transfer. The narrow hatched "instability 
strip" P p (M 2 i) bisecting the unhatched region is the plateau 
phase. This strip terminates at M 2 i ~ 3.5M , as for more 
massive donors M > Merit even when the radius expansion 
slows. 

We now consider the location of evolutionary sequences 
in Fig. 7, and discuss the relevance of P p and P cl -it for 
them. After turn-on of mass transfer a binary evolves 
to smaller M 2 and longer P. In the case of conserva- 
tive evolution the evolutionary tracks follow curves with 
(M - M 2 ) 3 Mf P = const. (M is the total binary mass), 
see the tracks of sequences S1-S3. In binaries with larger 
initial donor mass (M 2 i ii 2.5M ) the transfer rate is super- 
Eddington and the evolution non-conservative. The primary 
accretes at M e dd, and surplus material leaves the system at 
a rate M — Medd, carrying some specific angular momen- 
tum (of order the BH's specific orbital angular momentum) 
determined by details of how and where the material is ac- 
celerated. In this case the mass transfer rate is slightly lower 
and the period increase significantly faster than for conser- 
vative evolution, i.e. the evolutionary tracks would have a 
slightly larger slope |dP/dAf 2 | (cf. the example in Fig. 8). 

A simple corollary of (R2) is that a system will reach the 
transient plateau phase when its evolutionary track crosses 
the instability strip in Fig. 7. This is because P p is a unique 
function of the donor mass, independent of the initial donor 
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Figure 7. Orbital period — secondary mass (P — M2) plane for 
BHXBs, showing the exclusion zone for transients (unhatched). 
A BH mass of 8Mq was assumed. Systems in the unhatched re- 
gion are always persistent; systems in the narrow hatched insta- 
bility strip are always transient. Systems which are born in the 
large hatched area are always transient, but systems which have 
evolved from the unhatched into the large hatched area will re- 
main persistent for some time before they, too, become transient. 
Evolutionary tracks for sequences S1-S3 are also shown (solid 
where persistent, dashed where transient). The heavy dot along 
each sequence marks the location of the transient phase in the 
Hertzsprung gap. The dash-dotted line is the binary period when 
the secondary fills its Roche lobe on the turn-off main-sequence. 

mass. In contrast, as the initial donor mass dictates the 
mass transfer rate outside the plateau phase, we can ex- 
pect that the critical period where the system changes from 
a persistently bright to a transient X-ray source is close to 
fcrit(A/2i). Then, as a consequence of the opposite slopes 
of evolutionary tracks and the curve P cr it(A/2i) in Fig. 7, a 
system remains persistently bright when its track leaves the 
unhatched region, and becomes transient only later when its 
period roughly equals -P C rit(A/2i). The sequences S1-S3 (solid 
where persistent, broken where transient; heavy dots mark 
the plateau phase) and actual integrations of (^) confirm 
this. Hence no system in the unhatched region is transient, it 
is a "transient exclusion zone" , bisected by the Hertzsprung 
gap instability strip. The hatched region, on the other hand, 
may contain both persistently bright and transient systems. 

Fig. 8 shows how the transient exclusion zone depends 
on the BH mass M\. As M\ increases the steep branch of 
PmtiM?) moves to larger A/2, while the flat branch moves 
to longer P. The BH mass does not affect the location of 
the instability strip. 

5 SUMMARY AND DISCUSSION 

In extension of the work by KKB, KKS and King et al. 1997 
we studied the evolution of BHXBs with intermediate-mass 



Figure 8. As Figure 7, but with the critical boundary be- 
tween transient and persistent systems for different BH masses 
(dashed: 5Mq; solid: IOMq; dotted: 15Mq). The cross indicates 
the observed system parameters of GRO J1655-40 (Orosz & Bai- 
lyn 1997). The solid curves represent evolutionary tracks assum- 
ing conservative evolution (full line) and evolution with constant 
black hole mass (mass leaving the system carries the black hole's 
specific orbital angular momentum; dashed). 



giant donor stars in the context of the disc instability model 
for soft X-ray transients. These represent a hitherto some- 
what neglected group of systems with very high mass trans- 
fer rates and donors spending a significant fraction of the X- 
ray active time in crossing the Hertzsprung gap. We find that 
although systems in the process of crossing the Hertzsprung 
gap are in general persistently bright with Eddington- or 
super-Eddington transfer rates, there is a narrow instabil- 
ity strip where transient behaviour is possible. This strip 
extends over a secondary mass range 2.0 ^ M2/M-Q 3.5 
and is roughly given by P/d = 4.3M2/M0 —6.2, but the pre- 
cise location is subject to uncertainties in the stellar input 
physics. We predict that BHXBs in the Hertzsprung gap are 
not transient if the donor is more massive than ~ 3.5Mq, 
and that neutron star LMXB transients in the Hertzsprung 
gap do not exist. The latter is a consequence of a number of 
factors: in neutron star LMXBs the upper (stability) limit 
for the donor mass is smaller, the mass transfer rate is larger 
as the systems are closer to mass transfer instability, and 
the critical rate M cr it is possibly lower by almost an order 
of magnitude (see KKS). 

In our considerations we relate the evolutionary state 
of a system to the stability of its accretion disc, hence its 
appearance in X-rays, by comparing the evolutionary mean 
mass transfer rate M to the critical transfer rate M cr it for 
disc instability. Both the identification of the actual, instan- 
taneous mass transfer rate with M and the normalisation of 
M cr it are subject to some uncertainties: 

If the actual mass transfer rate deviated from the cal- 
culated evolutionary mean the system could of course be 
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transient where it is predicted to be persistently bright, and 
vice versa. Although such deviations from the secular mean 
rate are quite common in CVs (e.g. Warner 1987) there 
is no observational evidence for a similar large scatter of 
the mass transfer rate at a given orbital period in LMXBs 
(cf. van Paradijs 1996; Chen, Shrader & Livio 1997). Also, 
the most promising theoretical explanation for the observed 
spread of M in CVs is that the systems undergo mass trans- 
fer limit cycles caused by weak irradiation of the secondary 
(Ritter, Zhang & Kolb 1995; King et al. 1996). In LMXBs 
this mechanism does not work as the irradiating flux on the 
secondary is too large (cf. King 1998). Hence it seems not 
unlikely that the actual transfer rate is close to the evolu- 
tionary mean. 

The theoretical value of the critical rate Merit depends 
on a number of not very well determined quantities. Most 
notably these are the hydrogen ionization temperature in 
the disc, the relative disc thickness, the accretion efficiency 
r\ and the albedo. A major justification for the adopted nor- 
malisation comes from the fact that it matches the observed 
location of the critical rate separating transient from persis- 
tent sources (van Paradijs 1996). The remaining uncertainty 
in the normalisation of M cr i t does not affect the main results 
of this study, i.e. the existence and location of the transient 
instability strip in the Hertzsprung gap. This is a result of 
the steep gradient dM jdP of the transfer rate in this phase. 

The transient X-ray source GRO J1655-40 is the first 
confirmed member of the class of intermediate-mass BHXBs 
considered in this paper. GRO J1655-40 is also close to the 
SXT instability strip and, given the uncertainties in both the 
theoretical and observed values, might actually be a system 
right in the instability strip, consistent with its transient na- 
ture (e.g. Harmon et al. 1995, Tavani et al. 1996, Levine et 
al. 1996). Implications of this interpretation are discussed 
elsewhere (Kolb et al. 1997). Particular problems arise from 
the apparent lack of persistently bright systems with sim- 
ilar parameters in the Hertzsprung gap. It was suggested 
that these are not seen in X-rays as their effective photo- 
sphere might radiate at much longer wavelengths. Alterna- 
tively, GRO J1655-40 could be one of the predicted systems 
in the Hertzsprung gap with a stable accretion disc. It is 
by no means clear if such systems do appear as persistently 
bright sources. Instabilities in the accretion flow might cause 
dramatic changes in the effective photosphere, causing vari- 
ability in a given waveband. Kolb et al. (1997) suggested 
such a mechanism also for GRS 1915+105, the other X-ray 
source with apparent superluminal motion. 

A second possible intermediate-mass BHXB is 4U 1543- 
47. Recently, Orosz et al. (1998) determined its orbital pe- 
riod as 1.123 d and found that the components are likely to 
be a ~ 7Mq BH and a 2.1 — 2.5Mq main-sequence donor. 
This places the system below the turn-off main-sequence pe- 
riod (Ptms) line in Fig. 7 (dash-dotted; note that Ptms as 
shown in the equivalent Fig. 2 of Kolb et al. 1997 is slightly 
too small due to a calibration error). Therefore 4U 1543-47 
is, unlike GRO J1655-40, in a phase of Case A mass transfer 
where the transfer rate is determined by the donor's nuclear 
timescale. With Mi = 2.3M© and replacing K in (^) by the 
nuclear expansion dlnP/dr. on the main sequence we esti- 
mate M2 — 4 x 1O _9 M0 yr _1 , well below the critical rate 
KKS find for transient behaviour in BHXBs (but slightly 
above the critical rate claimed for neutron star systems, cf. 



KKB). A more detailed account of Case A mass transfer 
BHXBs is in preparation. 
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